Spin Polarization of Photoelectrons from Topological Insulators 
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We show that the degree of spin polarization of photoelectrons from the surface states of topo- 
logical insulators is 100 % if fully-polarized light is used as in typical photoemission measurements, 
and hence can be significantly higher than that of the initial state. Further, the spin orientation 
of these photoelectrons in general can also be very different from that of the initial surface state 
and is controlled by the photon polarization; a rich set of predicted phenomena have recently been 
confirmed by spin- and angle-resolved photoemission experiments. 



Three-dimensional topological insulators (TIs) arc 
strong spin-orbit interaction materials characterized by 
a bulk electronic gap and metallic topological surface 
state (TSS) bands with linear energy dispersions [lrl3|- 
The predicted linear energy dispersion of the TSSs in 
TIs were first observed in angle-resolved photoemission 
spectroscopy (ARPES) measurements TIs are con- 
sidered to be a promising candidate for spintronic devices 
because of their spin- momentum locking [H-0] . Aspects 
of the spin distribution of the TSSs in TIs have been 
measured by spin-resolved ARPES experiments [8l-fl3j|. 

Because the spin orientation of the photoelectron in the 
specific conditions used in previous spin-resolved ARPES 
studies agreed with the expected picture of the spin dis- 
tribution of the TSS electrons in TIs [Fig. 1(a)], photoe- 
mission matrix element effects were neglected in analyz- 
ing the spin polarization of the photoelectrons. (On the 
other hand, matrix element effects have been used in an- 
alyzing the circular dichroism of the TSS electrons in a 

ti yfir&-) 

Here, we find that the degree of spin polarization of 
photo-ejected electrons as defined in typical measure- 
ments with fully-polarized light is significantly higher 
(can in principle be 100%) than that of the TSSs (~ 50% 
for Bi2Se3 and Bi2Te3 according to first-principles cal- 
culations [6j ) , explaining the origin of high values of the 
measured degree of polarization (75% and > 85% in Pan 



et al. [12[ and Jozwiak et al. 13j, respectively, from ex- 
periments on Bi2Se3). Moreover, using the symmetries of 
the TI surface, we find that electron-photon interactions 
can completely alter the spin orientation of the photo- 
ejected electrons relative to that of the initial state and 
that the spin orientation of these photoelectrons can be 
controlled via light polarization tuning. For linearly po- 
larized light, the detected spin orientation is significantly 
altered except for the specific case of the wavevector k of 
an initial state being parallel to the in-plane component 
of the light polarization e. For k orthogonal to e, the two 
spins are predicted to be antiparallel to each other. More- 
over, for in-plane circularly polarized light, the spins of 
the photoelectrons are oriented cither completely parallel 
or completely antiparallel to the surface normal depend- 
ing on the handedness of the circular polarization. 



The measured degree of spin polarization is defined 
through the relation 



{*} k 



(1) 



where If and l_i are the intensities for the electron spin 
being aligned and anti-aligned with t, respectively. The 
unit vector t = i max which maximizes Eq. (flj defines 
the direction of the electron spin polarization. Hsieh et 
al. [9( reported in-plane P max [i. e. , restricting t in Eq. ([1]) 
in the surface plane] to be 20% for photoelectrons from 
the TSSs of Bi 2 Se 3 , Souma et al. [ijj and Xu et al. (n| 
reported 60% for those of Bi2Te3, and Pan et al. 12 1 
and Jozwiak et al. [l3| reported 75% and > 85%, re- 
spectively, for those of Bi2Se3. On the other hand, the 
degree of in-plane spin polarization of the TSSs for both 
Bi 2 Se 3 and Bi 2 Tc 3 obtained from first-principles calcu- 
lations is ~ 50% [6J. It is puzzling that the measured 
spin polarization of the photoelectrons 12, [l3| can be 
significantly higher than that of the corresponding TSS 
obtained from theory |6j. Due to, e. g. , spin-independent 
background signals, the measured degree of spin polar- 
ization is expected to always be lower than that from 
calculation. 

We find that the degree of spin polarization of photo- 
electrons from the TSSs of a TI is 100%. To show this, 
we start from the general expression for photocurrent 
from a detector (that selects a specific three-dimensional 
wavevector and hence a specific energy) with the spin 
quantization axis aligned with t: 



{f\E f =E,+hv} 



f) (f\n int \i) 



(2) 



where \i) is the initial TSS, |/) the photo-excited state, 
Ei and Ef their respective energies, hv the photon en- 
ergy, and H lnt the light-matter interaction Hamiltonian 
connecting the two states. State t, R_d) is the detected 
state: (i) its spin part is the eigenstate of s ■ t with eigen- 
value +1 (s is the Pauli matrix vector for spin half) and 
(ii) its spatial part is localized at Rd, where the detec- 
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tor is, far away from the sample surface. (Note that R^> 
is on the trajectory of the final-state wavepacket.) We 
dropped the obvious prefactor 5(E — E.- L — hv) in front of 
the right hand side of Eq. ([2]) due to the energy-resolving 
detector which collects electrons with energy E. A for- 
mal derivation of Eq. @ is reserved for interested read- 
ers [lU . Here we discuss the physical meaning of Eq. ([2} • 
First, we note that the detector reads spin character of 
the photo-excited state at Rd; hence, the near-surface 
part of the wavefunction affects the measured spin only 
indirectly through the matrix element (f\H mt \i). Sec- 
ond, a summation of the transition amplitude over de- 
generate photo-excited states |/)'s is necessary because, 
even in princi ple, we cannot tell which |/) is involved in 
the detection [19(. 

If we denote 

|/')= ]T \f){f\H^\i), (3) 

{f\E f =E l +hv} 

then we may rewrite Eq. © as 

Jfoc \(t,R D \ f')\ 2 ■ (4) 

Since the measurement is performed at Rd, far away 
from the sample where |/)'s are eigenstates of the free- 
electron Hamiltonian, the Bloch periodic part of the 
state |/') can be regarded as a position-independent two- 
element spinor in evaluating 7 t - by Eq. ([¥]). Then, from 
basic spin physics, we always can find a vector t' satisfy- 
ing |/') being the eigenstate of s • t' with eigenvalue +1 
and the eigenstate of s • (— t') with eigenvalue —1. Obvi- 
ously Ip and I_p = for this particular orientation 
and we have P ma x = 100% from Eq. fl}, i. e. , the degree 
of spin polarization of photo-ejected electrons is always 
100% regardless of that of the initial electronic state. 

An important ingredient that led us to this result is 
that the initial TSS electronic state \i) is not degenerate. 
If it is, a hole will be left in one of those different de- 
generate TSSs \i)'s after photodetection. Therefore, we 
can distinguish in principle which initial TSS is involved 
in the measurement; the detection probability amplitude 
corresponding to each \i) should be first squared and then 
summed, and not the other way round i.e., the 

photocurrent If will be the sum of contributions coming 
from all the degenerate initial TSSs |«)'s. If this happens, 
-Pmax 7^ 100%. The simplest example is, for a normal 
spin-degenerate material, If and I_f will always be the 
same regardless of the choice t, making P max = 0. 

Now we apply this general consideration to the case of 
the TSSs of a TI. According to recent ab initio calcula- 
tions [|| , the averaged degree of spin polarization of the 
TSSs, 

Plf = |<V>(n,k)|s|V>(n,k))| , (5) 



where \t/j(n, k)) is the two-component spinor wavefunc- 
tion, is roughly 50%. Because \ip(n, k)) is not an eigen- 
state of a spin operator t ■ s for any t, the degree of spin 
polarization for the TSSs had to be defined as an aver- 
aged quantity. On the contrary, the degree of spin po- 
larization of the photoelectrons P max [Eq. (P)] is 100% 
and that is so even if we take into account the imagi- 
nary part of the self energies of the involved electronic 
states because the detector probes the spin character (of 
the photoexcited states) at R^> instead of taking its av- 
erage over the entire space. If the imaginary part of the 
initial-state self energy is considered, the condition that 
"the initial electronic state is not degenerate" is naturally 
modified to "there are no other initial electronic states 
in the energy window within the width of the imaginary 
part of the initial-state self energy." 

This result tells us that a direct comparison between 
-Pmax [Eq. CQ)] and P^SS [Eq. ©] is not meaningful, 
and solves the app arent puzzle that the former from ex- 
periments J12J, |l3( is higher than the latter from theory 
(~ 50 %) [6| . Although the predicted degree of spin po- 
larization of the photoelectrons from the TSSs is 100%, 
the measured value will always be lower than this due, 
e. g. , to spin-unpolarized background signals and finite 
resolution of the apparatus. 

Equation ([2]) and our discussion based on it are not 
confined to topological insulators and can provide a guid- 
ance for the interpretation of any spin-resolved photoc- 
mission experiment using fully-polarized light. 

So far, we have discussed the magnitude of the spin 
polarization, without resorting to the details of the sys- 
tem. Now, by obtaining the state |/') in Eq. ((3]) using the 
specific symmetries of the TI surface, we determine the 
orientation of the spin polarization. We adopt the com- 
monly used effective Hamiltonian P^j(k) for the Bloch 
periodic part \<p(n, k)) of the wavefunction of TSSs in a 
TI with (in-plane) Bloch wavevector k = k x x + k y y (z is 
along the surface normal) given by 

i?xi(k) = hvk (sin 8^ a x — cos , (6) 

where v is the band velocity, 9^ the angle between k and 
the +k x direction, and o~ x and o~ y are the Pauli matrices 
acting on a two-component wavefunction, the so-called 
pseudospins. In constructing the effective Hamiltonian in 

Eq. ©, the basis states defining the ( J and ( ^ col- 
umn vectors, corresponding to pseudospin up and down 
states, are constructed from the two degenerate states at 
k = 0. 

Here, we consider a class of materials having the sym- 
metry of the surface of Bi 2 Se3 or Bi 2 Te3; however, the 
development can straightforwardly be extended to other 
classes. The two states at k = which are used as basis, 
\(f>i) and |02), can then be uniquely fixed by using the 
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FIG. 1: (a) The spin orientation of the surface states in a topological insulator on an equi-energy contour in the upper band 
[n = +1 in Eqs. ((8]) and ©]. (b)-(d): Similar quantities as in (a) for the measured photoelectrons with light having linear 



(along the +x direction), left-handed circular, and right-handed circular polarizations, respectively, 
handedness is defined from the viewpoint of the light source.) The surface normal direction is z. 



(For the latter two, the 



symmetry properties: 



e|0i) 

M\<h. 



101 



e|60 

M |0 2 ) 
Cs \6o) 



\<h) , (7) 



where is the time-reversal operator, M the reflection 
operator, x — ¥ —x (x is along the TK direction), and C3 
the operator for 27r/3 rotation around the z axis. The 
eigenvalue and Bloch periodic eigenstate are 



E(n, k) = n h v k . 



(8) 



and 
|0(n,k)) 



1 
75 



1) 



nie l6k \<j> 2 )) 



(9) 

respectively, where n = ±1 is the band index. The pseu- 
dospin expectation value of the TSS [Eq. ©] is thus given 
by 



(a) = n (sin 6*k x — cos 6*k y) 



(10) 



It is known that the actual spin expectation value (s) is 
aligned with the pseudospin expectation value Q, i. e. , 



(s) cx n (sin 9k x — cos 9^ y) 



(11) 



The orientation of the spinor eigenstates in the upper 
band (n = +1) is shown in Fig. 1(a). 

In calculating photoemission matrix elements, we use 
the theoretical framework of Wang et al. |15| . At small 
k, it is assumed that the final photoemission states 
|/)'s are spin-degenerate because they are well within 
the spin-degenerate bulk band continuum. In matrix 
element calculations for small k, we will approximate 
the periodic part of the final state Bloch wavefunctions 

|(k, kx)) and <^(k, k±)\ by those with k = and 



2m, 



(hv — En), where kj_ is the surface normal 
component of the photoelectron wavevector, m e the elec- 
tron mass and the energy of two- fold degenerate TSSs 
at k = 0. We denote these two (k-independent) states by 



and 



j_ / , respectively, and use the same symmetry 

relations as in Eq. ([7]) to define these states. Also, we 
will neglect the momentum dependence of the velocity 
operator v(k) = e~ lr ' k ve lr k [see Eq. (fl4)) ]. that has to 
be used in calculating the optical transition matrix ele- 
ment between the periodic parts of the Bloch wavefunc- 
tions. This theoretical setup [lH was employed to find 
the energy- and momentum-dependent spin polarization 
of the TSS in Bi 2 Se 3 from time-of-flight ARPES measure- 
ment with laser of energy 6.2 eV. The spin polarization 
of TSS electrons thus obtained is in excellent agreement 
with theory 120(1 and other spin-resolved ARPES experi- 
ments [8|, llOHl2j . Therefore, our theoretical development 
and predictions below should be valid when low-energy 
light source is employed. 



and 



An important point to note here is that 

are the actual spin-up and spin-down states along z far 
away from the surface in vacuum where the measurement 
is performed, because we have imposed the symmetry 
constraints of the system in Eq. ([7]) 0, |2fj . Even though 
the real spin character of those two basis states at the 
surface can be very complicated, we can regard the cho- 



sen degenerate doublet 



and 



the actual spin-up 



and spin-down states, respectively, from a measurement 
point of view. 

Next, the microscopic Hamiltonian Hd of an electron 
with spin-orbit coupling [21| is given by 



H n = 



2m e 



V 



Amlc 2 



(Wxp)-s, (12) 



where V is the one-electron potential. Using Peierls sub- 
stitution p — > p -A, where A is the vector potential 



4 



of the electromagnetic wave, and the relation (A) = 
Hd (p — #-AJ — Hd(p), we obtain the electron-photon 
interaction Hamiltonian 



where 



fl»t(A) = --A-v, 



V = — + ~A — 1T~2 ( S X W ) 



(13) 



(14) 



is the velocity operator plj ]. 

We first consider the matrix elements 

v Sjl = (0{|v|^) , 



(15) 



where ,s =t or J, and i = 1 or 2 are the index of the 
photo-excite state and the pseudospin basis index of the 
TSS state, respectively. If we define 



v± = v x ±iv v 



(16) 



and use the symmetry of the system for both the initial 
and final states [Eq. ([7}], only the following four are non- 
zero among twelve possible combinations [15j: 



V+ \<f> 2 



v- \4>i 



where a and f3 are real constants which can be determined 
from, e. g. , first-principles calculations. 

Plugging Eqs. (JT5J) and into Eq. flTHJl, we obtain 
the interaction Hamiltonian matrix H™f(A) connecting 
the two basis functions of the TSSs, \<f>\) and \<j>2), to 

the spin-up and spin-down photoexcited states, 



rrint 



(A) 



2c 



(Ay 



A x a y ) + i 



2/3 



A* I 



and 



(18) 

where / is the 2x2 identity matrix. In this study, we 
neglect the last term that depends on the z component of 
the light polarization, i. e. , we set A z = 0, in order to see 
the new physics clearly. Because this term is proportional 
to the identity matrix, it alone does not contribute to 
alteration of the spin direction of a photoemitted electron 
from that of the initial state. 

First we consider the case where the light is linearly 
polarized. Then, 



W$(A) = — eA (sin 9 A a x - cos A a y ), 



(19) 



where 9a is the angle between the vector potential A and 
the +x direction. As we discussed before, the photocur- 
rent is given by Eq. ©, where in our case \i) = \4>(n, k)) 



[Eq. ©] and |/)'sare 



and 



Since the photocur- 



rent / t - is nothing but the squared projection of the state 
I/') in Eq. ([3]) to the detector state If, Rd), it is essential 
to know the spin orientation of |/'). Using Eq. 

can write |/') [Eq. in the basis of 



and 



1 
71 



1 



(20) 



Comparing Eqs. © with O and using Eqs. dTUJ> 
and (ITT]) , we find that the net effect of photoexcitation 
on the detected electron spin polarization direction de- 
fined by the direction t maximizing (7 t - — 7_ t ~) / (i| + 
in Eq. (fTJ) is a rotation in direction through a change 



(21) 



where A#A.k = — ^k is the angle between the in-plane 
light polarization and the Bloch wavevector. The spin 
expectation value of the photoelectron arriving at the 
detector (s)^ in units of h/2 is thus given by 

(s) / = n [sin (6» k + 2A6» A .k) x - cos (6> k + 2A6» A ,k) y] ■ 

. (22) 

(Note that the magnitude of the spin polarization is 
100% in agreement with the above, general discussion.) 

For the special case of #k = 9a., the spin orientation of 
the photoelectron [Eq. (|22[) ] is the same as that of the ini- 
tial TSS electron [Eq. (| 1 1 [) ] . However, in general, the two 
are different. Especially, for states whose Bloch wavevec- 
tor k is perpendicular to the in-plane component of the 
light polarization A [i. e. , |A#A,k| = tt/2 in Eq. (|2"Tj) ]. the 
initial and final spins are antiparallel to each other. 

For in-plane circularly polarized light, the vector po- 
tential is given by A = A(x ± iy)/y/2 with the + and 
— signs denoting left-handed and right-handed circular 
polarizations as defined from the viewpoint of the light 
source, respectively. Then, within the effective Hamilto- 
nian formalism [Eq. (j!8[) ]. 



#ti(A) 



2V2( 



eA(±ia x - a y ) . 



(23) 



Applying the same argument as before, the spin polariza- 
tion direction of the photoelectron measured by a spin- 
resolved detector is that of 



|^'(n,k)) 



LHC 



and 



|0'(n,k)) 



RHC 



(24) 



(25) 



for left-handed and right-handed circular polarizations, 
respectively. The spin polarization direction of the pho- 
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toelectrons are thus pointed along the parallel (+z) and 
antiparallel {—£) directions to the surface normal, for 
left-handed [Fig. 1(c)] and right-handed [Fig. 1(d)] circu- 
lar polarized light, respectively. 

The phenomenon of photo-induced spin rotation in a 
TI predicted here was not observed previously. The rea- 
son is that A#A,k in Eq. (l22j) was held fixed or was 
allowed to change little in conventional spin-resolved 
ARPES measurements by linearly polarized light as dif- 
ferent k states were probed. Recently, Lanzara and 
coworkers 22j have tuned A#A,k and confirmed the 



photo-induced spin rotations, as predicted here for both 
linearly [Eq. d52J)] and circularly [Eqs. (JM]) and fl2SJ|] po- 
larized lights. 

In conclusion, we have shown that in spin-resolved pho- 
toemission experiments, the measured degree of spin po- 
larization of photoelectrons with a specific energy is 100% 
under ideal condition regardless of that of the initial state 
if the initial state is not degenerate and fully-polarized 
light is used. We have used this general principle to ex- 
plain why the degree of spin polarization of the photo- 
electrons from the surface states of a topological insula- 
tor can be higher than that of the initial states. Using 
the specific symmetries of the system, we have further 
shown that the spin polarization direction of photoex- 
cited electrons from the topological surface states is in 
general very different from that of the initial states and 
is dictated by the light polarization. Our results provide 
a theoretical basis for manipulation of the spin polariza- 
tion of the photoelectrons from the topological surface 
states [H-ii]. 
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SUPPLEMENTAL MATERIAL 
1. Fermi's golden rule 

Consider an electronic system described by a single- 
particle Hamiltonian Hq with energy eigenvalues and cor- 
responding eigenstates E n and \n), respectively, where 
n is the state index. Suppose now that we have 
a time-dependent perturbation of the form H lnt (t) = 
H lnt (A) e~ 2vwt , where A is the vector potential and hv 
the photon energy, in addition to the original Hamilto- 
nian Hq. The total Hamiltonian is given by 



H(t) = H a + H in \A)e 



(26) 



We are interested in the evolution of an electron which 
is at \i) before turning on the time-dependent perturba- 
tion. Assuming that the perturbation is turned on at 
t = 0, the wavefunction satisfying normalization condi- 
tion ((il)(t)\ip(t)) — 1) to first order can be written as 

\m) = e- lE * t/R \i) + e~ lE,t/h I/) , (27) 
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where c/(i) is the probability amplitude of finding an 
electron in state |/) at a given time t. Obviously, 



c/(0) = 0. 



(28) 



Using Eqs. $Z$, ([27]). and ([28"]). we can solve the time- 
dependent Schrodinger equation and arrive at 



where 



AEf cnin f L 

-' f -^(/| 

2 



AE f = (E f - Et) - hv. 



(29) 



(30) 
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Equation ([29]) leads to 



2. Spin-resolved photoemission experiment 



. 2 . AE f 

Pf(t)^\cf(t)\ 2 -^ S -^h\(f\ H*« 

Using the following relation 

lim =*(*), 



a->oo 7r a,X 



we finally obtain 



lim P f (t)/t 

t— ¥00 



7T . f AEf 



V 2h 



\(f\H 



int 



2tt 



5(AE f )\(f\ H» 



which is Fermi's golden rule. 



(31) 



(32) 



(33) 



Now we find an expression for photocurrent / t ~ in 
spin-resolved photoemission experiment with a detector 
aligned such that it captures electrons in a state t, R_d) 
of which (i) the spin part is the eigenvector of s • t with 
eigenvalue +1 and (ii) the spatial part is localized at the 
detector (r = Rd) as in the main manuscript. The pho- 
tocurrent is then given by 



I { oc lim K^RdI^C*))! ft. 



(34) 



Note that in the previous section, we have used |/) in 
the place of t, Rr>) to derive Fermi's golden rule; re- 
placement of the final state is the only difference. 

Using Eqs. and ([29]). we can write (i, R D \ip(t)) as 



(i,-R D \m) = E E - iel 

E {f\E f =E} 



AE, 



AE 
~2 



f- (/| H int \i) e- lE ^ h (t,R D \f) 



sini 



AE 

2h 



AE 

2 {f\E f =E} 



E (i,R D \f)(f\H"- 



(35) 



where AE = (E — Ei) — hv. An important knowledge the initial state is localized at the crystal and does not 
we made use of in Eq. ([35]) is that (t, Rd|«) = because extend to the detector. Obviously, a product of the terms 

in Eq. (|3"5j) corresponding to different E values does not 

' contribute to the photocurrent calculated from Eq. (|34p . 

Plugging Eq. J35| into Eq. flM|) and using Eq. ([32]) again, 
we obtain the expression for the photocurrent: 



J { cx E S ( E ~ E i - H 



oc 



E (i,R D \f) (f\ H" 

U\E f =E} 

2 



E (i,n D \f)(f\H ini \i) 

{f\E f =Ei+hv} 



r 



(36) 



which is Eq. (J2| of the main manuscript. Equivalently, Eq. (|36[) can also be written as 
I 



h (X 



E - ^ - H (t^Dlf) (/I # ir 



(37) 



